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ABSTRACT 



We present self-similar solutions for advection-dominated accretion flows with radial viscous 
force in the presence of outflows from the accretion flow or infall. The axisymmetric flow is treated 
in variables integrated over polar sections and the effects of infall and outflows on the accretion 
flow are parametrised for possible configurations compatible with the self-similar solution. Wc 
investigate the resulting accretion flows for three different viscosity laws and derive upper limits 
on the viscosity parameter a. In addition, we find a natural connection to non-rotating and 
spherical accretion with turbulent viscosity, which is assumed to persist even without differential 
rotation. Positive Bernoulli numbers for advection-dominated accretion allow a fraction of the 
gas to be expelled in an outflow and the upper limit on the viscosity predicts that outflows are 
inevitable for equations of state close to an ideal gas. 

Subject headings: accretion — hydrodynamics — Galaxyxenter 



1. Introduction 

Advection-dominated accretion flows (ADAFs) 
have been invented to explain low-luminosity black 
hole candidates like Sgr A* (Narayan et al. 1998) 
in our Galactic Center. The low-luminosity of 
this model is achieved in an optically thin plasma, 
where most of the energy is stored in hot ions, 
while electrons as potential radiators are ineffi- 
ciently coupled to the heat source and remain rel- 
atively cold. The electrons become nonetheless 
mildly relativistic close to the central black hole 
and the inevitable synchrotron radiation is ob- 
served from most ADAF candidates. The presence 
of magnetic fields not far from energy equiparti- 
tion with the gas is indicative of their origin in 
MHD-instabilities (Balbus & Hawley 1991) lead- 
ing to turbulence in the accretion flow and sub- 
sequent generation of an effective viscosity. On 
larger scales magnetic fields are likely to be re- 
sponsible for the collimation of outflows from ac- 
cretion disks into jets, seen in the cores of M87 
(Reynolds et al. 1996) and NGC 4258 (Lasota 
et al. 1996; Herrnstein et al. 1998), which 
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are prototypical ADAF candidates. Furthermore 
the model can explain the accretion in some low- 
luminosity AGNs of elliptical galaxies (di Matteo 
et al. 1999) and in NGC 4258 (Gammie, Narayan 
& Blandford 1999). In a recent X-ray survey 
(Sambruna, Eracleous & Mushotzky 1999) a few 
more examples have been found for low-luminosity 
core in radio-loud AGNs, which are candidates for 
advection-dominated accretion flows in their cen- 
tral engines, suggesting that ADAFs can be found 
even in radio-loud AGNs and that jet formation is 
a common feature. 

Outflow models for ADAFs have been investi- 
gated by Blandford & Begelman (1999) and ap- 
plied to several candidates (di Matteo et al. 1999; 
Quataert & Narayan, 1999). On the theoreti- 
cal side Igumenshchev & Abramowicz (1999) and 
Stone, Pringle & Begelman (1999) have performed 
time-dependent 2D calculations of accretion flows, 
which in some cases resemble ADAFs for certain 
viscosity parameters a 0.1, but suggest the pro- 
duction of outflows for larger a (Igumenshchev 
& Abramowicz 1999). It is found that the (rr) 
stress tensor component, which was not included 
in the original description of vertically integrated 
models for accretions flows, is important in the 
cited calculations for ADAFs. The existence of 
self-similar solutions with a radial viscous force 
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has been shown in previous work (Narayan & Yi 
1995a) and discussed for 2D solutions with a sep- 
aration of variables. 

In this paper we describe advection-dominated 
accretion flows in polar-integrated variables in- 
cluding the radial viscous braking force, which ei- 
ther produce outflows or are formed by wind infall. 
The wind infall is assumed to consist of free falling 
low angular momentum gas, which is cold with re- 
spect to the already existing accretion flow. One 
possible source for this gas are stellar winds from 
massive stars in young clusters as in the center of 
our galaxy. In the following we will refer to this 
wind infall into the accretion flow as infall. In 
our treatment the main difference between infall 
and outflow is the sign of the mass infall rate and 
we talk about winds, if it is not necessary to dis- 
tinguish between infall and outflow. We restrict 
the discussion to an extension of the self-similar 
solutions given by Narayan & Yi (1994) for the 
Newtonian limit. 

In §2 we present the equations, which describe 
the accretion flow including the reaction to winds. 
The role of (rr) stresses and bulk viscosity is em- 
phasised. We discuss angular momentum trans- 
port and viscosity in §3 and specify the possible 
equations of state and the resulting energy equa- 
tion in §4. General features of self-similar solu- 
tions are presented in §5 and in §6 detailed so- 
lutions for the a-viscosity law are shown. Con- 
sequences of the alternative /3-viscosity are dis- 
cussed in §7 and ADAFs with an intermediate 
shear-limited viscosity law follow in §8. We com- 
pare the solutions and draw our conclusions in §9. 

2. Stationary accretion with infall or out- 
flows 

In a first step we have to establish the set 
of equations, which describe the accretion flow. 
Advection-dominated flows, in which we are inter- 
ested, are known to be quasi-spherical (Narayan 
& Yi 1995a) and therefore we will use spherical 
coordinates in our discussion. Consider a station- 
ary, axisymmetric and rotating flow with angu- 
lar velocity around a compact object of mass 
M. Instabilities in the flow, either hydrodynamic 
or magneto-hydrodynamic in origin, generate tur- 
bulence on small scales and lead to an effective 
viscosity much larger than the microscopic one. 




Fig. 1. — Schematic view of accretion flow and 
associated outflow. The relativistic region around 
the central black hole (dashed line) is excluded 
from consideration. The axisymmetric flow is ro- 
tating with angular velocity Q and has a large ac- 
cretion velocity |tt|. The flow has an opening angel 
•& and polar-averaged quantities are averaged be- 
tween —1? and $ about the equator. The accretion 
flow may loose mass to an outflow as shown or 
acquire mass with an infall rate T,w 

The effective turbulent viscosity v in the flow will 
redistribute specific angular momentum £(r, 8) = 
sin 2 6r 2 Q by a local viscous torque between neigh- 
bouring rings or shells. Short term evolution on 
scales smaller than the mean free path of eddies A 
in the turbulent flow, which is related to the vis- 
cosity v = WoddyA is not resolved in this descrip- 
tion. All quantities like density or accretion veloc- 
ity must be understood in a local time averaged 
sense with probably large short term variations. 
Besides the time average, we will discuss accretion 
flows also as polar-averaged, one-dimensional flows 
with only a radial coordinate r. The average is 
taken over sections of shells occupied by the flow. 
Polar motions can consequently not be traced in 
this treatment. We are left with the angular ve- 
locity £1 = f2(r) and the radial component of the 
velocity, which is the accretion velocity u — u(r). 
Mass is accreted by the central object with a rate 
M = ~2irruT,. In this discussion S is a suitable 
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polar integral of the density 



,7T/2 

= 2/ 

Jir/2- 



d6 r sin dp 



(1) 



so that H = r-d is the vertical thickness of the 
accretion flow and the integral in equation (1) is 
restricted to a region around the symmetry plane 
of the flow to make room for winds seen in Fig.l. 
Nonetheless any vertical motion has a radial and 
polar component and the later is considered as a 
source or sink of mass, momentum and energy. We 
regard these sources as external to our description 
of the flow and call it infall into the accretion flow 
or an outflow respectively. 

For a non-relativistic flow the conservation of 
mass implies that the change in mass flux at every 
radius is balanced by mass exchange with the wind 
in the stationary case 
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where E^ is the mass per surface area added to 
the flow from the infall per time at a given radius. 
According to equation (1) this happens at a height 
n?, but rapid mixing with the flow is assumed so 
that the equations averaged over polar section oc- 
cupied by the accretion flow remain valid. Angu- 
lar momentum is a conserved quantity in accre- 
tion flows without winds and is redistributed by 
torques generated by viscous stresses t r< p 
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Here £ is a mass weighted polar average of £(r, 9). 
The infall contributes its own angular momentum 
to the flow Y^w^w and exerts an additional torque, 
which depends on the relative angular momentum 
of infall and flow £w — 

The accretion velocity u follows from the im- 
balance of gravitational force, centrifugal barrier 
and radial pressure gradient. In addition to that 
we include the viscous braking force F v from the 
(rr)-component of the stress tensor and bulk vis- 
cosity 
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with the integrated pressure V = J d6P and the 
radial velocity of the wind %, which adds its mo- 
mentum to the flow. The change in specific radial 
momentum induced by winds is proportional to 
the velocity difference Uw — u. The viscous force 
(see Appendix A) 
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has a contribution proportional the (rr) compo- 
nent of the shear tensor of the flow in equation (6) 
and the compression of the gas in equation (7). We 
make the crude assumption that the corresponding 
viscosities are the same for the (r<j>) and (rr) com- 
ponents of the shear stress and also for the bulk 
viscosity. Our ignorance is cast into the parame- 
ters r\\ and rj, which measure the strength of the 
(rr)-shear and bulk viscosity relative to the famil- 
iar (r</>)-shear viscosity respectively. For isotropic 
turbulence 771 w 1 is expected, but no estimate on 
77 can be derived from this assumption. We will 
call 77 and 771 viscous force measures in the follow- 
ing. Finally an energy equation must be given 
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which describes the change of specific internal en- 
ergy e of the flow due to compression (first term 
right side), viscous heating Q + , radiative cool- 
ing A, the enthalpy difference between flow w and 
wind oj\y and the kinetic energy associated with 
the velocity difference, which has to be dissipated 
into heat. Here v is the velocity vector of the flow 
and v-w the corresponding vector for the wind. 

3. Angular Momentum Transport and 
Viscosity 

The (rcf>) component of the stress tensor is as- 
sumed to be proportional to the corresponding 
component of the shear tensor with the kinematic 
viscosity fi = up as the factor of proportionality. 
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For the viscosity v we will adapt three different 
representations. The first and obvious one is the 
a-prescription introduced by Shakura & Sunyaev 
(1973) and the resulting a-disks will be discussed 
in §6. The a- viscosity is not a unique choice and 
effects of the so-called /3-viscosity introduced by 
Duschl, Strittmatter & Biermann are discussed as 
/3-disks in §(7). For the a- viscosity 
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K 



(9) 



we avoided to introduce the vertical scale-height 
H of the disk, because we believe that it leads 
to the misunderstanding that angular momentum 
transport in the radial direction depends on the 
thickness of the disk. Only with the assumption 
of vertical hydrostatic equilibrium is it possible to 
introduce H = c s /VLk in (9). 

The a-viscosity can be recovered from the (3- 
viscosity law in the case of shock-limited turbu- 
lence in a Keplerian disk. The /3-viscosity assumes 
that a typical length scale in the direction of trans- 
port Ar and the typical velocity difference Av^ 
between interacting shells, which exchange eddies, 
equal the mean free path of eddies and their typi- 
cal velocities relative to the mean flow respectively. 
A parametrisation with a new constant j3 is sug- 
gested 

v = Av^Ar = fiv^r . (10) 



If the typical velocity of eddies in a differen- 
tially rotating disk is limited by the sound speed, 
than the typical length scale Ar of communication 
between differentially rotating rings is estimated 
Duschl et al. to be 



(11) 



Only in disks with Keplerian rotation does this 
length scale equal the vertical scale-height and the 
form v — ac s H is recovered. As a third possibil- 
ity we consider the shear-limited form of the j3- 
viscosity with the length scale from equation (11), 
which has not been used for sub-keplerian accre- 
tion flows before. This is reasonable, if the mean 
free path of eddies is not limited by the vertical 
scale-height, but instead determined from the dis- 
tance between shells, for which the velocity differ- 
ence due to differential rotation equals the eddie 
velocity. This argument implies a larger effective 
viscosity for disks in sub-keplerian rotation. The 



(r<j>) component of the viscous stress tensor be- 
comes in any case 

dfl 



dr 



(12) 



We assume that the same viscosity prescription 
can be applied to the (rr)-component of the stress 
and the bulk viscosity. Their relative strength is 
scaled to the (r0)-viscosity by the force measures 
rji and r] introduced above. The heat generated 
by the described viscous forces has to be included 
in the energy equation and amounts to 

Q* = „x(,f' 2 ■ ' " 2 
\ or 



dr ( r ) 



+ 



d(r 2 u) 
dr 
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The internal energy is increased by viscous heat 
and compression and reduced by radiative cooling 
A. 

4. Equation of State 

In general accretion flows are expected to ex- 
ist for all reasonable equations of state for the ac- 
creted gas. The self-similar solutions we will use in 
what follows, are restricted to rather simple equa- 
tions of state. For ADAFs the gas should be well 
described by an ideal gas law with 7 = 5/3 for ra- 
tio of specific heats, because the temperatures in 
ADAF solutions are so large (T > 10 7 K for most 
radii of interest) that the gas is completely ionised 
and only ions of a few heavy elements retain their 
highly bound inner electrons (see Narayan & Ray- 
mond (1999) for possible emission lines of these 
elements). The change of the equation of state due 
to partial ionization can be neglected, but the con- 
tribution of magnetic to the total pressure must 
be considered for two reasons. The first is the 
observational evidence that at least some candi- 
dates (e.g. Sgr A*, see Narayan et al. (1998)) 
for ADAFs show a significant contribution of syn- 
chrotron emission to the total luminosity, which 
requires magnetic fields close to pressure equipar- 
tition with the gas. Given their existence, we 
might try a MHD-description of the accretion flow 
to separate the evolution of magnetic fields and 
gas, which might be different (Bisnovatyi-Kogan & 
Lovelace 1997) unless magnetic diffusivity (Hey- 
vaerts, Priest & Bardou 1996) and/or reconnec- 
tion require an even more complex model. The 
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easy way out is to assume that turbulence gener- 
ates small scale magnetic fields, which dominate 
the energy density in magnetic fields and produce 
an isotropic contribution to total pressure and en- 
ergy. In doing so, one arrives at a hydrodynamic 
description of the accretion flow with a equation of 
state, which has to incorporate the magnetic pres- 
sure. Consider an equation of state of the quite 
general form 

P ~ p Xp T XT ( 14 ) 

with the exponents XpiXT defined according to 
Cox & Giuli (1968). For the self-similar solutions 
we have to restrict the internal energy per unit 
volume to be proportional to the total pressure 



*q 7 -i)e 

XT 



(r 3 - l)e 



(15) 



with 7 = cp/cv being the ratio of specific heats 
and r 3 the adiabatic coefficient defined by Chan- 
drasekhar (1939). The restriction from equation 
(15) is that XpiXt and T 3 have to be constants, 
which requires that the relative strength of the 
magnetic pressure must be constant so that the 
effective equation of state is an ideal gas law with 
a ratio of specific heats less than 5/3. We define 
the isothermal sound speed c s = y^rV /T, from 
the integrated pressure and density. The sound 
speed used in the viscosity prescription is the to- 
tal pressure divided by the mass density. We use 
thermodynamic relations (Cox & Giuli 1968) to 
get an equation for the temperature from equation 
(8) for stationary accretion flows 

8T x P rc 2 s d (v\ , Q+-A 

+ -=r[ww-w + o ( 16 ) 



Here c v is the specific heat at constant volume 

(17) 



cv = 



X,(7-1)T 



and if we assume that gas pressure contributes a 
fraction 8 to the total pressure and magnetic fields 
are responsible for the rest, we get 

(18) 



\6 + (l-6) 



for the adiabatic exponent T 3 which is related to 
7 by equation (15). For constant Xp an d Xt we 



can rewrite the temperature gradient in equation 
(16) as a gradient of the sound speed and search 
for solutions in terms of the sound speed. 

5. Self-similar ADAF solutions with a 
wind 

The set of equations described in §2 and §3 al- 
low self-similar power-law solutions for accretion 
velocity, angular velocity and sound speed in the 
way Narayan & Yi (1994) have shown 



u = —UqS 



-v 2 c n = n aS ~^ 2 — (19) 
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with the radial coordinate scaled to the gravita- 
tional radius of the central mass 
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and the speed of light c. It is required that the 
cooling is either completely negligible or is a radius 
independent fraction of the heating rate 



/ = 1- 



A 
0+ 



(22) 



so that 1 — ./is the cooling efficiency, which will 
be small, if electrons are inefficiently coupled to 
the heat source. This happens, if ions are prefer- 
entially heated by viscous friction and the energy 
transfer rate between electrons and ions is smaller 
than the heating rate. 

For the self-similar solution the radial compo- 
nent of the wind velocity must be proportional 
to the accretion velocity and we get a measure 
£i = uw/u, which tells us the relative radial ve- 
locity of the wind, but we will neglect it in the fol- 
lowing discussion and assume £i = 0. This is mo- 
tivated from the infall calculations of Coker, Melia 
& Falcke (1999), which predicts a very small radial 
velocity of the infalling material in the disk mid- 
plane for the infall from stellar winds in the Galac- 
tic Center. If a fraction / larger than a few per 
cent of the viscously generated heat is not radiated 
away but kept in the flow as internal energy, the 
disk scale height H = c s /£Ik is of the same order 
as the radius and the infall may have a significant 
radial velocity component. Nonetheless we will 
neglect it in the following discussion. The same is 
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true for the enthalpy of the infalling material, but 
here we have better reasons to believe that the gas 
joining the ADAF is cold compared with the gas 
in the ADAF, which is heated by viscous friction, 
while the infalling gas maybe provided by stellar 
winds or the gas of the Hn region Sgr A West in 
case of the Galactic Center with temperatures of 
6 10 3 - 10 4 K (see Mezger, Duschl & Zylka (1996) 
for a review of the Galactic Center) . This should 
be compared to typical temperatures of ADAFs 
of 10 7 K at 10 4 black hole radii. The situation 
with outflows might be different, where it is hard 
to imagine that the gas leaving the disk has a dif- 
ferent temperature than the gas left in the flow. 
The gas might need some outward pointing radial 
momentum to get away from the ADAF, but again 
we will ignore it here. For infalling material we ex- 
pect it to be in free fall so that their total velocity 
is vw = vs = \/2cs _1 / 2 . In the case of outflows 
we require that they are able to escape the gravi- 
tational attraction of the central object and their 
velocity must be at least the free fall velocity vg. 

In the same way the rotation of the infall has to 
be a constant fraction £ = Q,y/ /Q of the rotation 
of the accretion flow itself and again calculations 
by Coker et al. (1999) for the Galactic Center 
predict infall with only weak rotation. Outflows 
on the one hand might be expected to carry their 
initial angular momentum from the point of origin 
and therefore £ = 1. On the other hand magnet- 
ically driven outflows will exert a torque on the 
remaining accretion flow if the gas in the outflow 
has different angular velocity than the point in the 
flow, to which it is connected by magnetic field 
lines. In this case < £ < 1 would be expected. 

The rate of mass added by the infall per surface 
area of the flow in its central plane is constrained 
to steep radial profiles and introduces one free pa- 
rameter p (Blandford & Begelman, 1999) not to 
be confused with the pressure P 

±w = -P— , S = S oS - 1 /2+p . (23) 
r 

For infall into accretion flows p has to be negative 
and the radial dependence of the infall rate must 
be steeper than r~ 2 . Therefore most of the mass 
is added at small radii and the total mass infall 
dM\y <~ dr r~ 1+p diverges in the center. This 
solution is consequently not valid for small radii, 
where the infall must deviate from the solution 



(23). Positive values of p correspond to outflows 
generated by the accretion flow and reversing the 
argument for the total mass now taken away from 
the ADAF, the solution is confined to small radii 
so that My/ is still smaller than the mass accre- 
tion rate supplied at the outer radius. The surface 
density of the flow implies a density p ~ /■ _3 / 2 +P ) 
which restricts the specific equation of state (14) 
to those with XpiXt and T 3 being constants and 
the temperature as a function of radius follows 

T XT _ r -l + (l-Xp)(-3/2+p) ^4) 

so that the radial pressure gradient is indepen- 
dent of the thermodynamic exponents and P ~ 
r -5/2+ P g ag density m S elf-similar solu- 

tions for advection-dominated accretion is uncon- 
strained and only the required inefficient cooling 
restricts the solutions to small mass accretion rates 
below w lCr 2 M E dd (Narayanfe Yi 1995b). M Edd 
is the Eddington accretion rate for an radiation 
efficiency of 10% in terms of the rest mass of the 
accreted material. 

6. a- ADAFs 

Whether the parametrisation of the wind, 
which is necessary for the self-similar solution in 
§5, is justified or not, depends on boundary con- 
ditions, which the wind has to meet at large dis- 
tances, or the internal physics of outflows. While 
outflows might naturally follow the self-similarity 
of ADAFs, the distribution of stellar wind sources 
will determine, if a power-law dependence of the 
infall rate is reasonable. The parametrisation of 
the viscosity (Shakura & Sunyaev 1973) in equa- 
tion (9) is based on dimensional arguments and 
besides numerical simulations, which do not exist 
for ADAF conditions, there is no way to deter- 
mine what values for a are appropriate in the case 
of stationary ADAFs. Without information on 
the strength of the magnetic fields from MHD- 
calculations, the adiabatic exponent T3 is also 
undetermined. We will discuss the self-similar so- 
lutions as functions of a and T 3 as well as for 
certain wind parameters and viscous force mea- 
sures T] and 771 . 

The self-similar ansatz in §5 reduces the dy- 
namical equations (3), (4) and (16) to a set of 
non-linear algebraic equations for the coefficients 
of sound speed a, angular velocity O and accre- 
tion velocity uq. The accretion velocity is derived 
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Fig. 2. — Isothermal sound speed c s for ADAFs in 
units of the local free fall velocity vg = y / 2GM/r. 
Both velocities have the same radial scaling and 
the ratio is independent of radius. Solutions are 
shown for T 3 = 1.4 as a function of the viscosity 
parameter a as long as the corresponding fl 2 is 
positive. The different sound speeds in all three 
cases of infall, no wind and outflow depend on the 
strength of the radial viscous force. The solutions 
correspond to values of rj and r\\ given in Figure 3 
for accretion and rotation velocity. 

from equation (3) in terms of sound speed a and 
a 

with a constant A, which reflects the torque ex- 
erted by the the wind on the flow. The combined 
wind parameter A is I for a vanishing wind or 
a non-rotating wind. For rotating infall the so- 
lution is restricted to moderate infall rates with 
p > —1/2. Otherwise accretion would not be pos- 
sible. For a given sound speed and accretion ve- 
locity, the centrifugal barrier and therefore the an- 
gular velocity follows from the radial momentum 
equation (4), if we know the viscous force mea- 
sures. We find for the square of the angular veloc- 



ity 3 

nS = c?-(|-p)o a -|a a ^ . (27) 

We are left with a quadratic equation for the 
square of the sound speed. One solution turns 
out to be irrelevant either by predicting a 2 to be 
negative or leaving Qq in equation(27) negative. 
But even the second root not always gives rea- 
sonable solutions, because with increasing a and 
small changes in a 2 as seen in Fig. 2 and Fig. 3, f2§ 
decreases and becomes negative. Solutions for the 
accretion problem exist only for values of the vis- 
cosity parameter smaller than a critical a c . The 
wind parameters for the solutions in Fig. 2 and 3 
have been chosen in a way that the infall or out- 
flow is in free fall or leaves with the escape speed. 
No radial momentum of the wind is included and 
both infall and outflow are cold. While the infall 
has no proper angular momentum, the outflow ro- 
tates with the SI of the accretion flow. With the 
requirements of a cold outflow and total velocity 
of the wind being the escape speed, which both 
enter the energy equation (16) in the same way, 
we make a minimal energy assumption for the ex- 
traction of internal energy by the outflow. For the 
equation of state we take the natural choice x p = 1 
so that the mix of gas and magnetic field behaves 
like an ideal gas. The ratio of specific heats equals 
T 3 , if we assume \t = 1- XT is not a parameter of 
our solutions, but determines the actual tempera- 
ture of the gas through equation (24). For energy 
equipartition we have T 3 = 1.4 from equation (18) 
and we use this value if not stated otherwise. For 
the cooling efficiency 1 — / we assume a low rate of 
1%. For the dynamics of the flow not the energy 
radiated away matters, but the energy left in the 
flow. The actual cooling efficiency is unimportant 
as long as it is small and / < 1. 

For an ADAF without winds it is possible to 
derive an analytic solution for the sound speed and 

3 

T = A{\- 2p(l - 6)) + r?(5 - 2p) + |t7i(1 + 2p) (26) 
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Fig. 3. — Radial accretion velocity — u in units of a 
times the local free fall velocity (left) and angular 
rotation velocity O of the flow relative to the lo- 
cal Keplerian rotation (right). Solutions are given 
for T3 = 1.4 and different strength of the radial 
viscous force. 77 measures the strength of bulk vis- 
cosity and rji the (rr) shear stress. Also shown is 
the standard ADAF solution without winds and 
no radial viscous force r\ = r]\ =0. The actual 
accretion velocity is almost linear in a and the 
difference at a — > is a difference in slope as a 
function of a. The rotation velocity Q, is shown 
only for ADAFs without a wind. 
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-1 -0.5 0.5 



p 

Fig. 4. — Excluded, allowed and marginally al- 
lowed regions of ADAF solutions with infall or out- 
flows in the plane of wind parameter p and ratio of 
specific heats, which equals T3 for \p = XT = 1 as 
mentioned in the text. O 2 becomes negative in the 
forbidden region and we can distinguish regions, 
where the Bernoulli number is either positive, rep- 
resenting unbound flows, or negative. Figure (a) 
shows the situation for corotation of wind and flow 
£ = 1. In figure (b) a non rotating wind £ = is 
assumed. The dividing line of allowed and forbid- 
den regions depend on the strength of the radial 
viscous force, for which 77 = and rji — 1 is as- 
sumed. For increasing a the allowed region moves 
deeper into the outflow regime and dividing lines 
are shown for a = 0.01, 0.1, 0.3, 0.6, 0.9. 
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determine the maximum allowed a c given by 4 
V8/9-2(r 3 -l)(7/3-r 3 )+B 



/(r 3 -l)(3r ? + 4r ?1 ) 



■ (29) 



In that case A equals 1 and the accretion velocity 
is uq = (3/2) a a 2 . If the viscosity parameter is 
small, the limit a — ► provides a good approxi- 
mation for the sound speed 



6/(r 3 - 1) 



3(r 3 -l)(5/-2 Xp ) + 10-6 Xp 



(30) 



and the angular velocity = \J c 2 — (5/2)a 2 as 
seen in Fig. 2 and Fig. 3. From equation (29) and 
(28) we immediately see that 1 < T 3 < 5/3 is re- 
quired for the no wind case and the familiar result 
arises that no solution exists for a non-relativistic 
ideal gas. It is also obvious that the upper limit 
on a increases and becomes irrelevant with van- 
ishing radial viscous force. In other words, with 
increasing viscous force measures r], rji ADAF so- 
lutions are restricted to reasonably small values of 
a < a c . 

One additional criterion for the realization of 
accretion solutions is the Bernoulli number 



GM 2 2 P 

Be = h r 2 Q 2 + u +e-\ 

r p 



(31) 



discussed by Narayan & Yi (1994, 1995a) for 
ADAFs and Blandford & Begelman (1999) for 
ADIOS. While the total energy decreases inwards 
the Bernoulli number being the total specific en- 
ergy plus P/p is positive and increases inwards 
for ADAFs. This changes with outflows due to 
their cooling effect by removing internal energy 
and reducing the Bernoulli number of the remain- 
ing flow. Correspondingly an infall increases the 
Bernoulli number of the flow. The combinations 
of allowed values for a and T 3 as a function of the 
wind strength p is shown in Fig. 4 for corotating 
(£ = 1) and non-rotating winds. The statement 



2 ( 5r? + - m 



(r 3 - 1) 



f*/Cr-i)(§-r, 



that ADAFs exist for an ideal gas with negative 
Bernoulli numbers, if an outflow is present (Bland- 
ford & Begelman 1999), is confirmed, even when 
the radial viscous force is included. The minimal 
possible wind strength depends on the viscosity 
parameter a as seen in Fig. 4. 

7. (3- ADAFs with winds 

Following the arguments in Duschl ct al. 
(2000) the viscosity law should not depend on 
the sound speed, if the actual velocity of turbu- 
lent eddies is smaller than the sound speed in the 
gas. In the /3-viscosity law (10) the eddy velocity 
is determined from the actual rotation velocity 
Woddy ~ y/Prft with the lessening yf/3 equally dis- 
tributed between velocity and length scale. The 
arguments given in §3 suggest that a-viscosity can 
be seen as the shock-limited case of the /3- viscosity 
and the a- ADAFs discussed in §6 can be checked 
for consistency, if this hypothesis is correct. Spe- 
cialising on the no wind case in the limit of small 
a we find the ratio 



c 2 s /M 
(rn)2 (5/3 - r 3Xp ) 



(32) 



(28) 



and the /3-viscosity can be applied for (3 < 1.5/ 
with the adaption of x p = 1 and T 3 = 1.4. There- 
fore the /3-viscosity is valid for ADAFs with low 
cooling efficiencies, if the resulting flows are simi- 
lar to a-disks. The scale length of the a-viscosity 
is the vertical scale height H = c s /VLk 1 which 
turns out be H w 1.55r v ///(1.6 + 6/) w 0.6r for 
a- ADAFs in the no wind case and the limit a — > 0. 
The length scale of the /3-viscosity is y/(3r and we 
expect the same sound speeds and angular veloc- 
ities for a- and /?- ADAFs below /3 < 0.3. This 
is verified by comparing Fig. 2 and 3 with Fig. 5. 
The actual solution for (3- ADAFs is derived in the 
same way as for a-ADAFs in §6 and we find the 
accretion velocity depending explicitly on and 
not c 2 . 

u = -^(3Q A (33) 

The constant A is defined in equation (25). The 
increased accretion velocity seen in Fig. 6 com- 
pared to a- ADAFs is explained by the larger value 
of Oo m equation (25) relative to the sound speed 
squared in case of the a-viscosity. The angular 
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velocity of the accretion flow can be given 5 

1 - (5/2 -p)a 2 
9/8[(l-2p(l-a)M 2 +W] [3 2 + l 



n 2 



(34) 



including the reaction of the flow due to the pres- 
ence of winds. Turning to the energy equation 
shows that two obvious solutions exist for the 
sound speed. First an non-rotating and non- 
accreting solution 



2c 2 



5-2p 



(35) 



where the gravitational attraction is completely 
balanced by the pressure gradient. The gas is in 
hydrostatic equilibrium and constitutes an opti- 
cally thin atmosphere. This is not a viable so- 
lution for a black hole, because general relativity 
does not allow a pressure supported atmosphere 
close to the horizon. The second solution is an ac- 
creting and rotating flow as shown in Fig. 5. The 
gas is cooler than in the atmosphere solution and 
strong outflows in corotation with the flow or even 
non-rotating outflows reduce the sound speed fur- 
ther so that the accreted gas can be gravitation- 
ally bound to the central mass. For a corotating, 
cold outflow, which reaches infinity with the es- 
cape velocity from the radius of origin £ 2 = 2, and 
a simple equation of state \p = 1 we can gi ye the 
sound speed for a case of negligible radial viscosity 
77 = = 0: 6 



/- [3(3/4 -p)f3 2 A 2 + l] P A 
(5/2 -p)/- 5 



(36) 



For non-rotating and cold infall, which is in free- 
fall for p < 0, or outflows barely reaching infinity 
(£ 2 = 1) with p < 0.36 the ratio of specific heats 
has to be smaller than a limiting value 



1 < 2 



l-p 



3 - 9p + 2p 2 



(37) 



H = ((5 - 2p)r t + 4/3(1 + 2p) m ) A 



^[4(6-p)(l-2p)(l + / 3 2 .4 2 ) 

/3 2 ^ 2 (3(l + 4p) 2 -p(17 + 2p))] 
.4(1 -p) (2 . 3 



+ 



r 3 



-/3M 2 (l-2p) 



0.45 



0.4 



0.35 



0.3 



0.25 



0.2 



0.2 



-,— , i 
p = -0.4 
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p - =.2:6. .No Rotation 



P = 0.6 



0.1 0.2 0.3 0.4 




p = -0.4 



0.1 0.2 0.3 0.4 



Fig. 5. — Isothermal sound speed relative to the 
local free fall velocity (left) and angular velocity 
in units of the local Keplerian value (right) for f3- 
ADAFs. For the no wind case p = the radial 
force measures arc rj = rji = 1 for the thick solid 
line, r\ = 0, rji = 1 for the long-short dashed line 
and no radial viscous force 77 = 771 = for the 
dotted line. The infall solution p = —0.4 has no 
rotation of the infalling gas and for the outflow 
solution with p — 0.6, the difference between a 
corotating outflow (thick solid line) and a non- 
rotating outflow is highlighted. 
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Fig. 6. — Accretion velocity (left) and Bernoulli 
number (right) for different /3-ADAFs. The iden- 
tification of solutions is the same as in Fig. 5. 



at which the rotation of the flow stops and no solu- 
tions exist for T 3 larger than that. This limit for 
T 3 equals 5/3 for ADAFs without winds or out- 
flows and is the well known result that ADAFs do 
not exits for an ideal equation of state. The con- 
sistency of /?- ADAFs follow from the above pre- 
dictions based on a-flows. With infall into the 
flow the sound speed increases and the rotation 
velocity decreases. The /3-viscosity law is there- 
fore more likely to be applicable to ADAFs with 
infall than with outflows. 

8. ADAFs with shear-limited viscosity 
and winds 

The motivation for a shear limited viscosity law 
is the observation that the standard a-viscosity 
agrees with the general /3-viscosity law only for Ke- 
plerian accretion disks with a shock-limited eddy 
velocity. In the shear-limited viscosity 



v — a 



n 



(38) 



we assume that the eddy velocity is indeed shock- 
limited and the mean free path of eddies is de- 
termined by the radial distance of differentially 
rotating shells, which allows the exchange of ed- 
dies along radial paths in a comoving frame of one 
shell. If the eddy velocity is the sound speed, the 
maximal distance of shells follows from a Taylor 
expansion of ft in Av^ 



Ar < ^ 

~ 2n 



(39) 



The length scale over which shells can interact 
is larger for sub-keplerian rotation. This is the 
main difference of equation (38) to the original 
a-viscosity. The consequence for the accretion ve- 
locity 



u 



3. a* 



(40) 



is an explicit dependence on the inverse of flo con- 
trary to (3- ADAFs in equation (33). The com- 
bined wind parameter A is defined in equation 
(25). The square of the angular velocity follows 
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from a quadratic equation with two roots 7 



1 5-2p 2 



5-2p 



-a 4 *. 
(41) 

In the limit of vanishing viscosity parameter a the 
angular velocity has two obvious solutions 



a 




(5/2 -p)a 2 



(42) 

For small a the rotation for the second solution 
is proportional to a and positive. The existence 
of these slowly or non-rotating ADAFs depends 
on non- vanishing force measures 77, 771 7^ 0. The 
solution disappears without a radial viscous force. 
The viscosity in non-rotating ADAFs tends to a 
finite value at a = 



1 - (5/2-rj)a 2 
a 2 * 



(43) 



and so does the accretion velocity. If we inter- 
pret dc s /f2 as the viscous length scale, it does not 
diverge for small a. The sound speed also tends 
to a positive value as seen in the upper solution 
branches for the sound speed in Fig. 8 and accre- 
tion velocity in Fig. 7 corresponding to the lower 
branches for the angular velocity in Fig. 7. The 
energy equation is itself quadratic in a 2 for either 
choices of ri in equation (41). Only one solution 
satisfies the condition 



a 2 < 



2p 



(44) 



for real values of f2. We end up with two solution 
branches, which join at a critical a c . Beyond a c 
no physical solution exists for the accretion prob- 
lem. Besides the slowly rotating solution discussed 
above, we find more familiar ADAFs, where the 
accretion velocity is linear in a as seen in Fig. 7. 
For small viscosity parameters, the solution is in- 
dependent of the force measures 77,771. Real solu- 
tions for the square of the sound speed in the no 



-a 2 A[(l-2 P (l-£ 1 ))A 
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-771(1 + 2p) + 77(5 -2p) 



0.4 



p 



0.2 



1 

- Outflow 


' 1 





-7) = 0, ^, = 0.5 






" No wind 




1 /' . 


'Wind infall 




~— -v. 

1 N 






/ 

/ - - 




, , 1 , 





0.1 



0.2 



0.3 



0.4 



0.8 



0.6 



c: 



0.4 



0.2 



~T»ind infall 




JI=77, = 



0.3 



0.4 



Fig. 7. — Accretion velocity (left) and rotation 
(right) for ADAFs with the shear-limited viscosity 
law. The ratio of specific heats equals T 3 = 1.4 
and wind parameters are p = 0.2 for the corotating 
outflow and p = —0.1 for infall with no angular 
momentum of the infalling material. The solution 
without radial viscous force is plotted as a solid 
thin line available for all values of a. Solutions 
with different viscous force measures are shown for 
the no wind case p = with 77 = 0, 771 = 1 as the 
thick solid line, 77 = 1, 771 = 1 the thin dash-dotted 
line and 77 = 0, 771 = 0.5 the short dashed line 
as indicated. For each solution the branch with 
large radial velocity shows slower rotation and vice 
versa. 
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Fig. 8. — Isothermal sound speed and Bernoulli 
number for the two branch ADAF model with 
shear-limited viscosity law. The identification of 
solutions is by line type and is the same as in Fig. 
7. The branch of each solution with large sound 
speed corresponds also to large accretion veloci- 
ties and slow rotation. The Bernoulli number for 
outflow solutions is smaller for the rapidly rotat- 
ing branch with small sound speed compared to 
the slowly rotating branch. In all other cases the 
Bernoulli number of the slowly rotating solution is 
smaller. 



a < 



(5/3- Xp r 3 ) 



3/(r 3 -lV»7 + (4/3)7/1 



(45) 



The situation does not differ from the standard a- 
viscosity in so far as no upper limit for a exists in 
the limit of vanishing radial viscosity r\ = r\\ = 
in the radial momentum equation (4) as long as 
the ratio of specific heats is in the range 1 < T 3 < 
5/3. In the limit of vanishing viscosity a — > the 
rotating solution coincides with the unique ADAF 
solution without winds and radial viscosity. For 
the simplest equation of state Xp = 1 the sound 
speed in the limit a — > is 



2e/ 



2e 2 + 5e/ 



5/3 - r 3 

r 3 -i 



(46) 



The critical a c , which sets the upper limit for the 
existence of ADAFs with shear-limited viscosity, 
depends not only on the force measures 77, 771, but 
also on the presence of winds. The same pattern as 
for the standard a- viscosity appears, in which so- 
lutions with outflows are possible for a wider range 
of a than for infall solutions as seen in Fig. 7 and 
Fig. 8. It turns out that the slowly rotating branch 
of solutions produces smaller Bernoulli numbers 
than the fast rotators for infall and the no wind 
case p — 0. Only for outflow solutions is the fast 
rotator preferred with smaller Bernoulli numbers. 
But in all cases shown in Fig. 8, the Bernoulli num- 
ber is positive and no gravitationally bound flow 
is found. This changes if more energy is extracted 
by the wind and the fast rotating outflow solution 
shown in the figures has negative Bernoulli num- 
bers, if the outflow is cold and reaches a terminal 
velocity equal to the escape speed at its origin. 
The allowed combinations of a and T 3 regardless 
of Bernoulli numbers are shown in Fig. 9. The di- 
viding lines in the (a, r 3 )-planc are vertical with- 
out a radial viscous force and bend in the way 
drawn in the figure, when the (rr) component of 
the stress tensor is included. 

9. Conclusions 

We have investigated the importance of a ra- 
dial viscous braking force for advection-dominated 
accretion flows (ADAFs) in the presence of in- 
fall or outflows. Under the assumption of almost 
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Fig. 9. — Dividing lines between regions with and 
without solutions for ADAFs in the a-T^ plane of 
viscosity parameter and ratio of specific heats for 
Xp = XT = 1 and only (rr) shear stress considered 
for the radial viscous force r) = 0, t]i = 1. The 
regions of allowed and excluded combinations of 
(a, r 3 ) depend on the wind parameter p, which 
varies from p = —0.3 to p = 0.4 in steps of Ap = 
0.1. The dividing line moves to larger a and 
for increasing wind parameter p. The no wind case 
p = is emphasised as the thick solid line. 

isotropic turbulence, bulk viscosity and (rr) com- 
ponent of the viscous stress tensor provide an ef- 
ficient brake of the rapid radial accretion of hot 
gas in ADAFs and produce additional heat due to 
viscous friction. This opens a second channel for 
transfer of kinetic into internal energy and sup- 
ports in part the radial motion so that the shear 
due to differential rotation and the centrifugal bar- 
rier is reduced in ADAFs with radial viscous brak- 
ing. 

We derived self-similar solutions of ADAFs for 
three different viscosity laws. The standard a- 
viscosity produces solutions, which show more and 
more sub-keplerian rotation with increasing a. 
The solutions terminate at a critical a c , which de- 
pends strongly on the ratio of specific heats so 
that the greatest possible a tends to zero for a 
non-relativistic ideal gas. At the critical a c the 
rotation of the flow vanishes and a purely radial 
inflow appears, which is not only supported by a 



pressure gradient, but also by the viscous braking 
force, provided the flow is still turbulent in the 
absence of differential rotation. This limit is rem- 
iniscent of Bondi accretion (Bondi 1952) in the 
presence of effective turbulent viscosity. The tran- 
sition from ADAFs to Bondi accretion nonetheless 
affords a suspicious fine tuning in a. 

The /3-viscosity law of Duschl et al. (2000) 
based on geometrical arguments in the absence 
of shock-limited turbulence allow ADAF solutions 
for all reasonable values of (3. The estimates of 
(3 as inverse of the critical Reynolds number of 
the flow suggests small /3s, for which the solutions 
do not differ from a-ADAFs. No transition to a 
Bondi like flow is possible in this case and no up- 
per bound on /3 exists. 

We showed that it is possible to derive a shear- 
limited viscosity law from the /3-viscosity men- 
tioned above, for which the transition to non- 
rotating Bondi like accretion with turbulent vis- 
cosity occurs naturally. The now familiar ADAF 
solutions with sub-keplerian rotation at a — > 
join with a second branch of solutions at a max- 
imal sustainable a. Similar to the a- ADAFs no 
accretion flow with larger a are possible. The sec- 
ond solution branch is a hot, slowly rotating, and 
rapidly accreting solution even for small values of 
a. The transition to viscous Bondi accretion oc- 
curs from the slowly rotating to the non-rotating 
solution in the limit a — > 0, where f2 is linear in 
a. A finite viscous force is present in these flows 
even in the limit a = 0. In all cases the connec- 
tion of ADAFs to non-rotating flows depends on 
the presence of a radial viscous force. 

The observation that ADAFs generally pos- 
sess positive Bernoulli numbers lead to the idea 
(Narayan & Yi 1995a; Blandford & Begelman 
1999) that ADAFs are good candidates for the 
production of outflows. In that way the accretion 
flow loses energy to the outflow and the remaining 
material is left with negative Bernoulli numbers 
and gravitationally bound to the central accreting 
mass. We confirm that statement in the presence 
of a radial viscous force for a- and /3- ADAFs and 
show the back-reaction of outflows on ADAFs for 
different outflow characteristics. Most noticeably 
is the cooling effect and the increased accretion 
velocity of the remaining ADAF as angular mo- 
mentum and internal energy is carried away. Out- 
flows with a minimal energy assumption for the 
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extracted energy have to be fairly massive p > 0.35 
to lower the Bernoulli number and leave a bound 
flow in case of the shear-limited viscosity. It is 
much easier to get a bound flow, if the minimal 
energy assumption is violated and the terminal ve- 
locity of the outflow is the escape speed from the 
origin of the out-flowing material. 

If a natural choice of the viscosity parameter — 
either a or a — exists and provided the strength 
of the radial viscous force can be estimated from 
isotropic turbulence, then outflows are inevitable 
for certain equations of state with ratios for spe- 
cific heats close to 5/3. This conclusion is inde- 
pendent of arguments based on the positiveness of 
Bernoulli numbers for ADAFs. 

The scenario of thin disk evaporation in binary 
systems (Liu et al. 1999) or the transition from 
cooling flows to ADAFs in low-luminosity cores of 
elliptical galaxies (Quataert & Narayan, 1999b) 
can explain the existence of ADAFs in these sys- 
tems. The formation of an ADAF, which is the 
most promising model for the spectral energy dis- 
tribution of Sgr A* in the Galactic Center, can- 
not proceed in either way. We suggest that an 
ADAF in the Galactic Center forms out of stellar 
wind infall (Coker et al. 1999) and the trans- 
fer of kinetic energy of the infall into internal en- 
ergy of an advection-dominated flow. The even 
larger Bernoulli numbers produced in this way 
wound give rise to subsequent outflows and re- 
duce the mass accretion rate inferred from in- 
fall calculations by Coker et al. (1999) to the 
smaller accretion rates predicted from spectral fit- 
ting (Quataert & Narayan, 1999) of ADAF mod- 
els to Sgr A*. 

I thank Ramesh Narayan and Wolfgang Duschl 
for helpful conversations and challenging discus- 
sions. This work has been supported through 
DAAD fellowship D/98/27005. 
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A. The radial viscous force 



REFERENCES 



The radial viscous force in spherical coordinates 
is derived from the tensor divergence of the shear 
stress tensor tij and the bulk viscous force from 
the divergence of bulk viscosity and compression. 
The stress tensor is assumed to be proportional to 
the shear tensor V which is defined as 



1 



1 



{vi-j+Vj-i) - -gijV i 



(Al) 



with the kinematic viscosity \i = vp as a scalar 
function relating shear and stress tij — fiDij . Here 
Vi are velocity vector components, gij the metric 
tensor and Vi-j implies covariant differentiation of 
the velocity. The radial viscous force derived from 
the shear tensor is 



4 fd_ ' 
3 \dr 



(A2) 



d /u\ a /u\ 

dr \r J dr \r ) 



d 



where axial-symmetry and vg = is assumed and 
u is the radial velocity as used throughout the pa- 
per. The corresponding heating rate due to this 
viscous friction is 

2 



(<9r (r) 



(A3) 



The contribution from compression and bulk vis- 
cosity C to the radial viscous force is 



/bulk = (Cv.l);r = 



d f( d 



{ur 2 ) 



dr \r 2 dr 
and the corresponding heating rate is 



9bulk 



C d 2 
^d~r {Ur ] 



(A4) 



(A5) 



Multiplying by r and performing polar integration 
of the force gives the radial viscous force used in 
(6) and (7) 
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and 



d (vY, d . 2n 
F^ = r-{^-(ur ) 



(A6) 



(A7) 



with the replacement of the bulk viscosity by the 
effective turbulent viscosity v. The heating rates 
used in equation (13) are derived from equation 
(A3) and (A5) in the same way. The force mea- 
sures r], 7/i used in the text have been omitted here. 
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